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A 29 nucleotide deletion in open reading frame 8 (ORF8) is the most obvious genetic change in severe 
acute respiratory syndrome coronavirus (SARS-CoV) during its emergence in humans. In spite of 
intense study, it remains unclear whether the deletion actually reflects adaptation to humans. Here we 
engineered full, partially deleted (—29 nt), and fully deleted ORF8 into a SARS-CoV infectious cDNA 
clone, strain Frankfurt-1. Replication of the resulting viruses was compared in primate cell cultures as 
well as Rhinolophus bat cells made permissive for SARS-CoV replication by lentiviral transduction of 
the human angiotensin-converting enzyme 2 receptor. Cells from cotton rat, goat, and sheep provided 
control scenarios that represent host systems in which SARS-CoV is neither endemic nor epidemic. 
Independent of the cell system, the truncation of ORF8 (29 nt deletion) decreased replication up to 23- 
fold. The effect was independent of the type | interferon response. The 29 nt deletion in SARS-CoV is 

a deleterious mutation acquired along the initial human-to-human transmission chain. The resulting 
loss of fitness may be due to a founder effect, which has rarely been documented in processes of viral 
emergence. These results have important implications for the retrospective assessment of the threat 
posed by SARS. 


Emerging zoonotic viruses such as Severe Acute Respiratory Syndrome coronavirus (SARS-CoV) are a consider- 
able concern for public health’. Bats of the genus Rhinolophus are the bona-fide animal reservoir for SARS- and 
SARS-related CoV (SARSr-CoV)’*. Epidemics may emerge from wild animal reservoirs in which a plethora of 
viral variants exists. After initial cross-host infection, the occurrence of positive selection with adaptive changes 
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is thought to be essential for viral emergence**. The confirmation of phenotypic changes acquired in the process 
of adaptation requires virus isolation, which is not always possible. In the case of SARS-CoV, isolation of virus 
strains from the bat reservoir has widely failed’, except in singular instances involving bat-borne viruses whose 
spike protein showed an exceptionally high amino acid identity with the epidemic strain®*. Additional reasons 
for low isolation success include the lack of appropriate cell culture systems as well as low virus concentrations 
and often inappropriate storage conditions for samples from field investigations. Studies of phenotypic properties 
of reservoir-borne viruses can be facilitated by the reconstruction of genetic traits through reverse genetics. 

In the present study we reconstructed variant SARS-CoVs carrying different forms of open reading frame 
(ORF) 8, an accessory gene in the SARS-CoV genome that is among the most variable genes in bat-associated 
SARSr-CoV. ORFS is also one of the most relevant genes when studying potential viral adaptation to humans, 
as the ORF8 coding sequence has undergone gradual deletion during the human epidemic. Early epidemic 
SARS-CoVs contained a full ORF8 that is also present in the genomes of almost all bat- and carnivore-associated 
precursor viruses’"'*. A 29 nucleotide (nt) deletion within ORF8 occurred in all strains involved in the middle 
and late phase of the human epidemic”’’. Subtotal or total deletions of up to 415 nts occurred within and around 
the ORF8 region in SARS-CoV strains circulating in the very late phase of the epidemic. We have previously 
found that a European SARS-related CoV carried by rhinolophid bats in Bulgaria did not contain any ORF8 or 
similar gene” while virtually all SARS-related CoV from Asia possess a single, continuous ORF8*'*””. 

The 29 nt deletion in ORF8 was the most obvious genetic change during human-to-human transmission of 
SARS, causing the expression of truncated gene products termed ORF8a and ORFS8b. It has been widely hypoth- 
esized that the truncated products led to a modulation of pathogenicity or replication that favored adaptation 
of SARS-CoV to humans”'’””. For instance, it was found that replication of SARS-CoV is increased in cells that 
overexpress the protein encoded by ORF8a"°. In contrast, the same protein led to an attenuation of replication 
when engineered into recombinant infectious bronchitis virus (IBV)'’. The 8b protein, when overexpressed, 
was found to induce apoptosis and to be involved in cellular degradation of the viral envelope protein during 
SARS-CoV infection'®'®!’, When engineered into an IBV infectious clone, protein 8b inhibited the induction of 
interferon (IFN) as triggered by poly(I:C)-stimulation and thus enhanced IBV replication. This effect was found 
to involve binding of ubiquitin and interference with induction of the IFN type 1 promoter via IRF3”. Another 
study found that the full 122 amino acid protein encoded by ORF8 induces ATF6-dependent transcription, which 
triggers the expression of chaperones and leads to a general attenuation of the protein translation level, thus mod- 
ifying the unfolded protein response”. Loss of ORF8’s original cellular localization in the endoplasmatic retic- 
ulum would have ablated this function in human viruses of the middle and late epidemic phases”””'. However, 
the effect that this change had on SARS-CoV replication level remains unclear. Studies of SARS-CoV replication 
in experimental mice have not provided indications as to essential functions of ORF8a or 8b”. However, these 
studies involved high inoculated virus doses and utilized a mouse model that may not reflect relative changes 
of replication in transition from bats to humans. In sum, the available data leave it unclear whether the changes 
that occurred in ORF8 led to a modification of viral replication when comparing bat versus human hosts, and 
specifically, whether the deletion of 29 nt involved an increase of replication level in human cells as often impli- 
cated. Genome deletions have also been observed in MERS-CoV, and also here it was speculated that deletion 
may reflect adaptation to the human host or the release of selective pressure exerted exclusively in the zoonotic 
reservoir”, 

In the present study on SARS-CoV, we first followed up on our initial finding of the absence of ORF8 in 
genomes of SARSr-CoV from European bats, and sequenced rhinolophid bat-associated SARSr-CoV from four 
countries across Europe (Spain, Bulgaria, Italy, Slovenia). We then studied the influence of ORF8 on replication 
in the context of a full replicating SARS-CoV genome, based on a host transition scenario represented by cell cul- 
tures. We constructed recombinant viruses with full ORF8, truncated ORF8 (29 nt deletion), as well as completely 
deleted ORF8. Replication was compared in primate cell lines (VeroFM, MA104), a novel cell line generated from 
the lung of a rhinolophid bat, three additional non-chiropteran cell lines, as well as human airway epithelial cul- 
tures. We find that the 29 nt deletion conferred an attenuation of replication level irrespective of the host system 
studied. 


Results 

Bats in Europe carry SARS-related CoV that lack ORF8. We have previously described the detection of 
SARS-related CoV in Rhinolophus in Bulgaria'’. The full-length genome of the virus sequenced in that study con- 
tained no ORF8. To further investigate the presence of ORF8 in SARS-related CoV in Europe, we analyzed fecal 
samples from rhinolophid bats in Spain, Bulgaria, Italy, as well as Slovenia. All five Rhinolophus species commonly 
encountered in Europe were tested positive for coronavirus RNA (Table 1). However, none of 92 geographically 
and phylogenetically representative viruses provided evidence for the presence of an ORF8 gene. 19 of these 
viruses were fully sequenced, confirming absence of ORF8 at orthotopic or heterotopic genome positions. 


Generation of recombinant SARS-CoV encoding ORF8 variants. To understand the consequences 
of the absence of ORF8, recombinant SARS-CoV (rSCV) were generated as shown in Fig. la. The variant termed 
rSCV rg contained a complete ORF8 as encountered in bats, civet cats, as well as in early human cases of SARS, 
thus resembling the pre-epidemic and starting phase of the SARS outbreak. The variant termed rSCV,,; corre- 
sponds to isolates from the main phase (also referred to as middle and late phases) of the epidemic containing a 
29 nt deletion that introduces a ribosomal frameshift separating ORF8 into ORF8a and ORF8b. Two alternative 
virus variants without ORF8 were constructed. Variant 1, rSCVg.js3_;, had ORF8 removed without any replace- 
ment. Variant 2, rSCV4,s_>, had ORF8 replaced by a 5 nt spacer sequence (AATAA) which occurs instead of ORF8 
between ORF7b and the nucleocapsid gene in a natural SARS-related CoV variant from European rhinolophid 
bats’. 
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R. euryale (n = 44) 44 
Bulgaria ene n ae (n= 1) 

R. mehelyi (n= 4) 4 
Italy 45 R. ferrumequinum (n=6) 6 
Slovenia”? | 36 R. hipposideros (n=5) 5 
Spain 285 R. hipposideros (n=21) 21 
Total 872 All 5 species known in Europe | 92 (100%) 








Table 1. Bat specimens included in this study and confirmation of absence of ORF8. Absence of ORF8 
was determined by RT-PCR using reverse primers F29260R (TTTGTATGCGTCAATGTGCTTG) for 
reverse transcription and F28182R (GGGTCCACCAAATGTAATGCGG) and forward primer F27626F 
(GAGAAAGACAGAATGAATGAGC) for PCR. 
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Figure 1. Generation and evaluation of ORF8 variant recombinant SARS-CoV. Variants of the open reading 
frame 8 (ORF8) were designed in accordance to their appearance in nature. (a) rSCV¢¢,, represents a single 
ORF$8 as found in reservoir bats and amplification hosts in China, as well as in the early phase of the SARS 
epidemic. The middle phase of the epidemic was dominated by a virus variant carrying a 29 nt deletion 
resulting in the disruption of ORF8 in 2 reading frames, 8a and 8b, as seen in rSCV,,;. The absence of ORF$8 is a 
genomic feature of reservoir bats in Europe and the epidemic virus in the late phase. Two deletion variants were 
constructed. Variant 1, rSCV4,is_;, perfectly misses ORF8. In variant 2, rSCV4e83-2 ORF8 is replaced by the short 
substitutional sequence AATAA in accordance to the upstream region of the nucleocapsid gene of the European 
SARS-related bat-CoV BtCoV/BM48-31/Rhi bla/Bulgaria/2008 (NC_014470). (b) Plaque morphology 

of rSCV sfn and rSCV.,; were very similar, while rSCVaeg-1 produced only diffuse and rSCV4eas-2 reduced 
plaques. Western Blot analysis revealed that only rSCV stun rSCVep and rSCVaeg-2 infected cells expressed 

the nucleocapsid protein, while none could be detected in cells infected with rSCV4eag-1. Detection of 8 actin 
served as loading control. Virus replication of the three ORF8 variants, (rSCVgeig = rSC Vaeig_2) Was monitored 
by plaque titration after infection of VeroFM cells at two different multiplicities of infection, 1 (c) and 0.001 (d). 
Virus growth was determined in at least 3 independent experiments in triplicates. Shown is one representative 
experiment. Error bars represent standard deviation of the mean. 


Infectious viruses were rescued from cDNA as described”. 24h after co-electroporation of in-vitro transcribed 
RNA and nucleocapsid subgenomic RNA, parts of the cell culture supernatant were transferred to VeroFM 
cells and virus replication was monitored by real-time RT-PCR. 72h post infection (hpi), rSCV ge, rSCVepi and 
rSCVieg-2 gave evidence of replication by real-time RT-PCR, while rSCV4eg-1 RNA increased more slowly, reach- 
ing a plateau level only by 96 hpi. 
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Figure 2. Replication of ORF8 variants in primate cell cultures with and without IFN pre-treatment. (a) 
VeroFM and MA104 cells were infected with three ORF8 variant viruses at MOI = 0.001. Supernatants were 
harvested at 24 hpi and virus titers determined by plaque titration. The experiment was done in triplicates. 
Significance of replication differences between virus variants was determined by t- Test (***p < 0.001, 

*p < 0.05). (b) VeroFM cells were incubated with universal type I IFN at indicated concentrations 16h prior to 
infection with rSCV gfu and rSCV4eg at MOI=0.001. At 24 hpi supernatants were plaque titered. Experiments 
were done at least twice in triplicates. Shown is one representative experiment each. Error bars represent 
standard deviation of the mean. 
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All variants were plaque-quantified. During this process, differences in plaque morphology were noted. While 
rSCV gga and rSCV,,; yielded plaques of similar sizes, rSCV4eg-2 Showed smaller plaque size and rSCVaeg-1 gener- 
ated no distinct plaques (Fig. 1b). 

The absence of plaques in rSCV4eg-ı could point to a replication defect caused by the loss of essential gene func- 
tions, including the nucleocapsid gene whose expression is indispensable for virus replication*®. Modifications 
of ORF8 affect the upstream context of the transcription regulatory sequence (TRS) of the nucleocapsid gene”. 
Western blot analyses showed that rSCV ge, rSCVepi» as Well as rSCV4eis-2 expressed similar amounts of the nucle- 
ocapsid protein, while no expression was seen in cells infected with rSCV4eg-ı (Fig. 1b). Subsequent experiments 
were therefore carried out using variant 2, which is hereafter referred to as rSCV4,is. 


Integrity of ORF8 facilitates the replication of recombinant SARS-CoV. Replication of rSCV gfu 
rSCV,,;, and rSCVg.ig were tested in VeroFM cells at high and low multiplicity of infection (MOI) (1 and 0.001 
plaque forming unit (PFU)/cell). Supernatants were harvested at 8, 24, 48, and 72 hpi and titered in VeroFM cells 
(Fig. 1c, d). At MOI= 1, 30-40% of the cells were already lysed by 24 hpi, precluding meaningful comparisons 
of viruses (Fig. 1c). At MOI=0.001, replication was in the exponential phase by 24 hpi (Fig. 1d). rSCVg¢,1, grew 
significantly more efficiently than both other variants at 24 hpi (p < 0.05 and p < 0.001 comparing to rSCV.,; and 
rSCVaeig, respectively, using t-Test). The fully deleted variant rSCV4e replicated significantly less efficiently than 
rSCV,,; at 24 hpi (p = 0.002, t-Test). By 48 hpi, replication of all variants reached plateau levels. 

The replication facilitating effect of ORF8 is IFN independent. Due to the essential role of the type I 
IFN system in the restriction of virus infection, it was tested whether the ORF8-dependent replication phenotype 
might be linked to the type I IFN response. MA104 monkey kidney cells were used because they are competent 
for IFN induction and signaling”. VeroFM cells would only reflect differences in IFN signaling, due to their 
defect in IFN beta gene expression”. According to previous results, virus infections with all three variants were 
carried out at MOI=0.001 and virus progeny was measured at 24 hpi by plaque assay. 

All viruses generally replicated more efficiently in VeroFM cells than MA104 cells. However, relative differ- 
ences between strains were very similar in both cell lines, suggesting that IFN beta gene induction that is only 
functional in MA104 cells may not play a significant role for the observed phenotypes (Fig. 2a). 

To test the potential influence of the ORF8 gene product on IFN signaling, virus growth after external addition 
of IFN was compared. To this end, VeroFM cells were incubated with indicated concentrations of universal type 
I IFN alpha prior to infection with rSCV gry and rSCV gae. At 24h p.i., rSCV gr. grew to titers 30 times higher than 
rSCV4geig (Fig. 2b). If cells were treated with IFN prior to infection the overall replication of both viruses decreased 
in a dose-dependent manner. However, the replication difference between both viruses remained constant. This 
suggests that ORF8 has a SARS-CoV-specific replication-enhancing effect that is independent of IFN antagonism. 


ORF8 improves replication of recombinant SARS-CoV in a reservoir bat cell line. In previous studies, 
it has been suggested that the function of ORF8 may be more relevant for replication in cells from the actual 
animal reservoir, than for replication in human or primate cell cultures*®. Because bats of the genus Rhinolophus 
are the bona-fide animal reservoir for SARS-related CoV’, an epithelial cell line was generated from the lung 
of Rhinolophus alcyone (Fig. 3a). One pregnant female individual was euthanized and embryonal tissues were 
prepared using techniques described previously (Fig. 3b)°'~*’. The resulting immortalized Rhinolophus lung cell 
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Figure 3. Generation of an hACE2-transgenic bat cell line for infection experiments with ORF8 variant rSCV. 
(a) Picture of an African Rhinolophus alcyone bat (copyright Victor Max Corman). Primary cell culture and 
immortalization of R. alcyone embryonic lung cells (b) was done as described in the Material and Methods 
section. To determine the RhiLu cell type the epithelial protein marker cytokeratin (c) and the fibroblast marker 
S-100A4 (d) were detected by immunofluorescence assay using mouse-anti cytokeratin or rabbit-anti S-100A4 
Ig. Secondary detection was performed by incubating with goat-anti-mouse cyanin 2- or goat-anti-rabbit 
cyanin 3-labeled Igs. The bars represent 20 um. (e) Integration of hACE2 into the genome of RhiLu cells was 
verified by PCR. The vector containing the hACE2-puromycin resistance gene construct was used as a PCR 
positive control. (f) Expression of hACE2 was confirmed by Western blot analysis using mouse anti-hACE2 Ig 
(1:1,000). In addition, 8 actin protein was detected using rabbit anti-B-actin Ig (1:2,000) to ensure that similar 
protein amounts were applied. MA104 cells, expressing hACE2 naturally, served as positive control. (g) Virus 
replication of rSCV ge, rSCVepi and rSCVaeg was observed by titration of supernatants sampled at 24, 48, and 

72 hpi. Cells were infected in triplicates at an MOI of 0.001. Error bars represent standard deviation of the mean. 
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culture, termed RhiLu, clearly expressed cytokeratin as a marker for epithelium (Fig. 3c), and hardly showed any 
expression of the fibroblast marker S-100A4 (Fig. 3d), suggesting that the cell culture mainly contains epithelial 
cells. Initial infection experiments in RhiLu cells identified these cells to be non-permissive for SARS-CoV infec- 
tion. By use of VSV-G protein-pseudotyped lentiviral particles, RhiLu cells were therefore transduced with the 
gene for the SARS-CoV entry receptor, human angiotensin-converting enzyme 2 (hACE2), and selected to stably 
express hACE2 upon co-expression of a puromycin resistance gene from the same lentiviral vector. Chromosomal 
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Figure 4. Replication of ORF8 variants in different “non-SARS-CoV-host” cell lines and differentiated human 
airway epithelial cells. (a) Cotton rat, goat, and sheep cells were transduced with lentiviruses to transiently 
express the SARS-CoV receptor hACE2 for at least 72 h. Expression of hACE2 was verified by Western Blot 
analysis, detection of 8 actin served as a loading control. To improve clarity blots were cropped. Full length blots 
are presented in Supplementary Fig. S1. (b) Cells were infected in triplicates with different ORF8 variants at 
MOI 0.001 24h after lentiviral transduction, 24 hpi supernatants were sampled and virus replication determined 
by plaque titration. (c) Differentiated human airway epithelial cells were infected in triplicates at MOI=0.1, 
supernatants were sampled at 48 and 72 hpi for plaque titration. The experiment was done twice. Shown is one 
representative experiment. Error bars represent standard deviation of the mean. Significance of replication 
differences between virus variants was determined by t- Test (***p < 0.001, **p < 0.001, *p < 0.05, n.s., not 
significant p > 0.05). 
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integration of the hACE2 gene and expression of the hACE2 protein were verified by RT-PCR and Western blot 
analysis (Fig. 3e, f). The resulting hACE2-transgenic cells termed RhiLu-hACE2 were highly permissive for 
SARS-CoV infection. 

RhiLu-ACE2 cells were infected with all three ORF8 virus variants at MOI = 0.001. Supernatants were har- 
vested 24, 48 and 72 hpi and titered. As shown in Fig. 3g, the same hierarchy of viral replication levels as in 
VeroFM cells was observed (rSCV grun > rSCVepi > rSC Vacig). These results did not point to a specific effect of the 
gene product of ORF8 on replication in bat cells but rather suggested a general promoting effect on virus replica- 
tion conferred by ORF8. 


SARS-CoV ORF8 determines replication in non-host cell lines and in human airway epithelial 
cultures. To determine whether the effect was more general and independent of host taxa, a low-passage air- 
way epithelial cell line from Sigmodon hispidus (cotton rat)** and primary lung cell lines from Capra hircus (goat) 
and Ovis aries (sheep) were transduced using vesicular stromatitis virus G protein pseudotyped lentiviruses to 
transiently express hACE2 as described above (Fig. 4a), and subsequently infected with the ORF8 virus vari- 
ants at MOI = 0.001. Again the same hierarchy of viral replication levels, namely rSCV stun > rSCVepi > rSCVaeiss 
was observed at 24 hpi (Fig. 4b). Virus rSCV¢¢, replicated significantly less efficiently than rSCV,,; in all three 
“non-host” cell lines (t-Test: p = 0.001, p < 0.05, p < 0.01 for cotton rat, goat and sheep, respectively). Replication 
differences were levelled off by 48 hpi, when all viruses reached the plateau of replication. Of note, the utilized cell 
cultures in total represent hosts from four different orders of Euarchontoglires and Laurasiatheria, the two main 
groups of mammals. 

As the virus variants constructed in this study corresponded to viruses from the three major phases of the 
human SARS epidemic, we were interested to see whether the deletions have an effect on human respiratory tract 
infection as represented by in-vitro differentiated human airway epithelial cultures (HAE). Infection was done at 
MOI= 0.01, and virus production was measured by plaque titration after 48 and 72 hpi in supernatants. At both 
time points rSCV sfn replicated to significantly higher titers than rSCV,,,; or rSCV4eg in HAE (Fig. 4c). 


Discussion 
Here we have shown by viral reverse genetics and advanced cell culture models that SARS-CoV ORF8 facilitates 
viral replication irrespective of the host cell system. These data suggest the occurrence of an attenuating mutation 
in the initial phase of the human SARS epidemic. 

SARSr-CoV is a paradigmatic pathogen for the study of viral reservoirs and the processes involved in epi- 
demic emergence®*'*'>°. The acquisition of certain spike proteins by recombination may have formed the 
viral lineage that emerged from the reservoir and established itself in humans during the SARS epidemic®®. 
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The gradual deletion of ORF8 constituted the most obvious change in SARS-CoV after emergence. Because few 
changes occurred in other parts of the genome, we have focused our study on ORF8 and kept the rest of the 
genome constant representing the sequence of SARS-CoV strain Frankfurt-1. This prototype virus was isolated 
from the late phase of the SARS-CoV epidemic when transmission chains already occurred in countries outside 
China’”®. Even in this late epidemic strain, a reconstitution of the full ORF8 reading frame led to an increased 
replicative capability, suggesting that the viral genome was still compatible with its primordial ORF8 element in 
spite of possible onward evolution in humans’. 

Compared to earlier observations with ORF8-deleted SARS-CoV, our results show a clear phenotypic differ- 
ence in replication in relevant models of human respiratory tract infection. The phenotypic difference depends 
on the presence of full ORF8, with enhanced replication as opposed to the deletion variant. Earlier experiments 
in which no clear differences in replication have been noted were aimed at the identification of strong attenuation 
markers as necessary for live vaccine development, and hence worked at very high MOIs at or above 1”. We have 
also included such high MOIs in our experiments for comparison, and our results were highly concordant with 
those earlier studies. For instance, like in the study by Yount et al.”*, differences between viruses with full ORF8 
and 29nt deletions were only about three-fold at MOI = 1. Natural infection, however, does not involve high virus 
doses*’. When a virus is passed from human to human, the stochastic nature of infection success implicates that 
inocula just below or above one unit of human-infectious virus are transmitted. For instance, during influenza A 
transmission in ferrets as few as 2 virus units were transmitted between animals**. Under conditions of inter-host 
transmission, the observed phenotypic differences as observed in our study may have significant effects on viral 
fitness”. 

In endemic viral infections, loss of fitness should cause viral lineage extinction in competition with more 
reproductively capable lineages that co-circulate within the host population’’. After a single-time zoonotic intro- 
duction, however, competing viral lineages are unlikely to exist. Variants with slightly deleterious mutations, ran- 
domly selected through transmission bottlenecks, can continue to reproduce in spite of reduced fitness - an effect 
known as founder effect. Interestingly, it has been shown by in-vitro studies that virus populations with reduced 
initial fitness suffer less from slightly deleterious mutations than populations that replicate on peak fitness level“. 
Reduced initial fitness is a condition that can be expected in early-stage zoonotic epidemics when the virus is not 
yet adapted to the new host environment’. 

Considering the present results, we therefore suggest that the 29 nt deletion in SARS-CoV is the result of 
a founder effect that has permitted survival in spite of reduction of fitness. This interpretation contrasts with 
the earlier notion that the 29 nt deletion reflects adaptation to humans. The conclusion of adaptation is partly 
based on results from expression of ORF8 or its truncation products in overexpression systems or expression in 
heterologous virus genomes, suggesting various influences on virus-cell interaction’®’. For instance, one study 
provided evidence for IFN evasion mediated by ORF8b"’. This is not confirmed by our experiments studying 
ORFS in full virus context. Other authors have proposed that the 29 nt deletion may have been neutral for fit- 
ness after viral host transition to humans, assuming that ORF8 may elicit a function that is only relevant in the 
bat host”’*. However, we have not observed any bat-specific effects in cell culture, and rather show that ORF8 
optimizes fitness irrespective of the host cell system, including hosts that are irrelevant for the SARS-CoV chain 
of emergence’’. Our experiments suggest that only the deletion of further portions after initial fragmentation of 
ORF8 may have been neutral to fitness. Such deletions were seen during very late phases of the SARS epidemic 
in Hong Kong’. 

Deletions in accessory reading frames were also oberserved in MERS-CoV’*******. Transmission of deleted 
variants was confirmed in an outbreak in Jordan, involving deletions in ORF4a, ORF3 and potentially other 
parts of the genome’. The available studies leave it open whether these deletions involved changes of repli- 
cation level or virulence. However, it is known that ORF4a acts as an effective antagonist of MDA5-dependent 
induction of type I IFN“, and that MERS-CoV is highly sensitive against type I IFN in human airway epithelial 
cultures”. Based on these known mechanisms, attenuation rather than human adaptation should be considered. 
Our results are also relevant in the context of a series of experimental studies recently conducted to understand 
potential adaptive changes in the 2014 Ebola virus Makona outbreak in West Africa. Whereas initial studies 
suggested human adaptation with increase of replication level during the outbreak*°®, later studies found that 
a late-outbreak strain rather caused reduced virulence and prolonged survival of experimental animals”’. These 
results provide another reminder of the fact that outbreak-associated mutations do not have to increase replica- 
tion or virulence. 

It is interesting to consider the consequences for viral propagation conferred by a reduction of viral replication 
level. As pointed out in Marzi et al. for Ebola virus Makona, a prolonged survival time as seen with late epidemic 
strains could have increased the duration of infectious virus shedding in humans and may thus have increased 
the long-term fitness of those viruses”. We cannot, at present, exclude whether similar effects may have pro- 
vided a fitness advantage to SARS-CoV on host population level, such as by keeping infected individuals socially 
interactive for prolonged times while infected with a slightly attenuated virus. On an individual level, however, 
a reduction of replication level is likely to attenuate the pathogenicity of infection and reduce the health burden 
caused by a given outbreak. 

It may be seen as a weakness in our present study that no experimental animals were infected as in Marzi et al. 
However, mouse models of SARS-CoV do not reflect human disease as accurate as macaques do for Ebola virus, 
and the HAE culture system used in the present study already provides an appropriate model for the authentic 
site of replication of SARS-CoV in the human body (Ebola virus infection cannot be modeled by organ-specific 
cultures). Moreover, the infection phenotypes as seen in our study for ORF8,,, and ORF8,,; have already been 
demonstrated in mice, and the corresponding effects in cell culture based on those variants have been reproduced 
in our study”. We have therefore avoided additional animal experimentation and focused on human epithelial 
models with low inoculation doses such as seen in natural infections. 
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Our data suggest that SARS-CoV has suffered an attenuating mutation by the 29 nt deletion that constitutes a 
landmark genetic change. The SARS epidemic in 2003 may have taken a more severe course if not involving this 
mutation. Further work, including work in experimental animals, will be required to understand and confirm 
whether the absence of ORF8 in European bat-associated SARSr-CoV correctly predicts lesser epidemic risks as 
compared to Asian strains. 


Materials and Methods 

Sample collection and processing. In total, 827 bats were sampled in four countries (Bulgaria, Italy, 
Slovenia, Spain). All animals were handled according to national and European legislation for the protection 
of animals (EU council directive 86/609/EEC). Licenses for sampling of bats using mist nets, hand nets or harp 
traps were obtained from the respective countries and authorities: Bulgarian Ministry of Environment and 
Water, permit No. 192/26.03.2009"°, Italien Ministry of the Environment, permit No. 192/26.03.2009°, Slovenian 
Environment Agency, permit No. 35701-80/2004°°, Service for the Biodiversity Conservation of the Rural 
Counseling of the Xunta de Galicia, Spain, permit No. 52/2010 n.s. 13697. No animals were sacrificed during this 
study. All animal handling and sampling was done by trained personnel, with animal safety and comfort as the 
first priority during minimally invasive sampling (collection of faeces). Bat species were identified on site and, if 
necessary, mitochondrial DNA in representative fecal samples was amplified and sequenced for species confirma- 
tion as described previously”. Captured bats were freed from nets immediately and put into cotton bags for 2 to 
15 min to allow them to calm down before examination. While being kept in bags, bats produced fecal pellets that 
were transferred to 500 ul RNAlater RNA stabilization solution (Qiagen, Hilden, Germany) for sample processing. 
After homogenization, 50 ul of the suspension was resuspended in 560 ul of buffer AVL from the Qiagen viral 
RNA minikit and processed according to manufacturers instructions. The elution volume was 50 ul. 


General cell culture procedures. Cells were grown in Dulbecco’s Modified Eagles Medium (DMEM) as 
described earlier*!. For titration of rsCV Vero E6 cells (ATCC CRL-1586) were used. For infection studies African 
green monkey kidney cells, VeroFM (kindly provided by Jindrich Cinatl, University of Frankfurt) and MA104 
(a gift from Friedemann Weber, University of Marburg), Rhinolophus alcyone (R.alcyone) embryonic lung cells 
(RhiLu, prepared in-house as described below), as well as bronchial epithelial cells from Sigmodon hispidus (cot- 
ton rat) and lung cells from Capra hircus (domestic goat) and Ovis aries (domestic sheep) were used””®. HAE 
cultures were generated and cultured as described elsewhere”. 


Generation and characterization of R. alcyone lung cell cultures. Bats were caught in Ghana under 
research permit no. CHRPE49/09; A04957 Wildlife Division, Forestry Commission, Accra, Ghana. Primary bat 
cell culture and immortalization of RhiLu cells by lentiviral transduction of the simian virus 40 large T antigen and 
genotyping were done as previously described*'*’. To determine the RhiLu cell type the epithelial protein marker 
cytokeratin and the fibroblast marker S-100A4 (calcium binding protein A4 or fibroblast specific protein 1) 
were stained by immunofluorescence assay using mouse anti-cytokeratin (ab7753) or rabbit anti-S-100A4 immu- 
noglobulins (Ig, ab27957; both supplied by abcam, Cambridge, UK). Secondary detection was performed by 
incubation with goat anti-mouse cyanin 2- or goat anti-rabbit cyanin 3-labeled Igs. 


Generation of a SARS-CoV susceptible bat cell line by lentiviral transduction. Since RhiLu cells were 
not susceptible to SARS-CoV the receptor hACE2 was stably transfected by lentiviral transduction”. Genomic 
integration of hACE was verified by PCR on genomic DNA. PCR was performed using a hACE2 gene specific for- 
ward primer (GAATGTAAGGCCACTGCTCAACTA) and a puromycin resistance gene specific reverse primer 
(TCAGGCACCGGGCTTGC) yielding a 1.9kb amplicon. Expression of hACE2 was confirmed by Western blot 
analysis as described earlier**”’. Protein lysates were separated on a 10% SDS-PAGE gel, and blotted onto a 
0.45 um polyvinylidene fluoride membrane. Primary detection of hACE2 was done using a mouse anti-hACE2 Ig 
(1:1,000; R&D Systems, Wiesbaden-Nordenstadt, Germany), for secondary detection a goat anti-mouse horse- 
radish peroxidase (HRP)-conjugated Ig (1:20,000) and SuperSignal® West Femto Chemiluminescence Substrate 
(Fisher Scientific, Schwerte, Germany) were used. As a loading control samples were analyzed for 8 actin expres- 
sion with a rabbit anti-Bactin Ig (1:2,000; Sigma-Aldrich, Munich, Germany) and a goat anti-rabbit horseradish 
peroxidase-labeled Ig (1:20,000)°". 


Generation of recombinant SARS-CoV. Recombinant SCV was generated as previously described”. 
The single full-length ORF8 was generated by inserting 29 nts into ORF8a/b. In an overlap extension PCR two 
templates, generated by using primers F26020F (CGGCTCTTCAGGAGTTGCTA) in combination with 29nt-rev 
(ITCCATTCAGGTTGGTAACCAGTAGGACAAGGATCTTCAAGCACATGA) and 29nt-fwd (CTGGTTA 
CCAACCTGAATGGAATATAAGGTACAACACTAGGGGTAATACT) with pB-fwd (GCCCTTAAACGCCT 
GGTTGCTAC), were fused. Underlined nts indicate overlapping regions leading to the introduction of 
29 nts. The PCR product was inserted into subclone pEF via restriction sites BamHI and NotI. ORF8 was deleted 
from subclone pDEF by Phusion® Site-Directed Mutagenesis (Fisher Scientific) using 5’-phosphorylated 
primers delO8-1-fwd (AATGTCTGATAATGGACCCCAATCAAACCAACGTAGTGC) and delO8-1-rev 
(GTTCGTTTAGACTTTGGTACAAGGTTCTTCTAGATCC) or delO8-2-fwd (LAAAATGTCTGATAATGGA 
CCCCAATCAAACCAACG) and delO8-2-rev (TTGTTCGTTTAGACTTTGGTACAAGGTTCTTCTAGATCC). 
Underlined nts are the substitutional sequence for ORF8 compared to delO8-1. Assembly of full-length 
SARS-CoV genome plasmids and rescue of recombinant viruses were done as described before’. Briefly, 
full-length SARS-CoV plasmid was linearized by NotI and in-vitro transcribed (mMESSAGE mMACHINE® 
Kit, Applied Biosystems, Darmstadt, Germany). Capped mRNA was electroporated into baby hamster kidney 
cells and supernatant was subsequently transferred to susceptible VeroFM cell culture 24h post electropora- 
tion. Recombinant virus was harvested three days post infection. The ORF8 mutations were verified by PCR 
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using a specific reverse primer F29260R (TTTGTATGCGTCAATGTGCTTG) for reverse transcription and 
ORF8 covering forward primer F27626F (GAGAAAGACAGAATGAATGAGC) and reverse primer F28182R 
(GGGTCCACCAAATGTAATGCGG) for conventional PCR. Sequencing of the PCR product ensured the 
integrity of the introduced mutations. Expression of the nucleocpasid was confirmed by Western blot analysis as 
described above. Primary detection was done using a rabbit anti-nucleocapsid Ig (1: 500; abcam). For secondary 
detection a goat anti-rabbit HRP-conjugated Ig (1:20,000) and and SuperSignal® West Femto Chemiluminescence 
Substrate (Fisher Scientific) was used. 


Generation of transiently hACE2 expressing cell lines. Cells were seeded according to their size and 
growth rate to yield 80% confluence in 24-well plates. After attachment overnight, cells were infected with lenti- 
viruses to yield 50 ng reverse transcriptase activity per well in a reduced cultivation volume of 200 uL DMEM per 
well for at least 24h. Expression levels of hACE2 were determined by Western blot analysis as described above in 
a time course experiment. Protein expression levels were constant for three days after transduction. Therefore, 
cells were infected with rSCV 24h post transduction as described below. 


Virus infection. Cells (4 x 10° cells/mL) were seeded in 24-well plates and if necessary pre-incubated with 
200 ul recombinant universal type I IFN alpha (PBL InterferonSource, Piscataway, USA) for 16h prior to infec- 
tion. Infections with rSCV gftn rSCV,,; and rSCVq.ig were done at an MOI ranging from 0.001 to 1. Virus was 
diluted in OptiPRO™ serum-free medium or in HBSS buffer in case of HAE cultures. Cells were inoculated for 
1 to 2h at 37°C, washed twice with PBS or HBSS, and supplied with fresh DMEM or HAE cell culture medium”. 
Supernatants were taken at designated time points usually 8, 24, 48, and 72h post infection and stored at —70°C 
for titration or real-time RT-PCR analysis. All virus containing samples were mixed with equal volumes of 0.5% 
gelatin in OptiPRO™ (stock solution 5% gelatin in water) for stabilization of infectious particles. All infec- 
tion experiments were done under biosafetly level 3 conditions with enhanced respiratory personal protection 
equipment. 


Plaque titration. Titration of rSCV was done as previously described”. Vero E6 cells (3.5 x 10° cells/mL) 
were infected with a serial dilution (in OptiPRO™) of virus infected cell culture supernatants for 1h at 37°C. 
After removing the inoculum cells were overlaid with 2.4% Avicel (FMC BioPolymers, Brussels, Belgium) 1:2 
diluted in 2 x DMEM supplemented with 2% Penicillin/Streptomycin, 2% L-glutamine, 2% non-essential amino 
acids, 2% sodium pyruvate and 20% fetal bovine serum. Three days after infection the overlay was discarded, 
cells were fixed in 6% formaldehyde and stained with a 0.2% crystal violet, 2% ethanol and 10% formaldehyde 
containing solution. 

The datasets generated during and/or analysed during the current study are available from the corresponding 
author on reasonable request. 


References 
1. Jones, K. E. et al. Global trends in emerging infectious diseases. Nature 451, 990-993 (2008). 
2. Drexler, J. F., Corman, V. M. & Drosten, C. Ecology, evolution and classification of bat coronaviruses in the aftermath of SARS. 
Antiviral research 101, 45-56 (2014). 
3. Wolfe, N. D., Dunavan, C. P. & Diamond, J. Origins of major human infectious diseases. Nature 447, 279-283 (2007). 
4. Holmes, E. C., Dudas, G., Rambaut, A. & Andersen, K. G. The evolution of Ebola virus: Insights from the 2013-2016 epidemic. 
Nature 538, 193-200 (2016). 
5. Graham, R. L., Donaldson, E. F. & Baric, R. S. A decade after SARS: strategies for controlling emerging coronaviruses. Nature 
reviews. Microbiology 11, 836-848 (2013). 
6. Ge, X. Y. et al. Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor. Nature 503, 535-538 
(2013). 
7. Yang, X. L. et al. Isolation and Characterization of a Novel Bat Coronavirus Closely Related to the Direct Progenitor of Severe Acute 
Respiratory Syndrome Coronavirus. Journal of virology 90, 3253-3256 (2015). 
8. Hu, B. et al. Discovery of a rich gene pool of bat SARS-related coronaviruses provides new insights into the origin of SARS 
coronavirus. PLoS pathogens 13, e1006698 (2017). 
9. Molecular evolution of the SARS coronavirus during the course of the SARS epidemic in China. Science 303, 1666-1669 (2004). 
10. Chiu, R. W. et al. Tracing SARS-coronavirus variant with large genomic deletion. Emerging infectious diseases 11, 168-170 (2005). 
11. Guan, Y. et al. Isolation and characterization of viruses related to the SARS coronavirus from animals in southern China. Science 
302, 276-278 (2003). 

12. Lau, S. K. et al. Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats. Proceedings of the National 
Academy of Sciences of the United States of America 102, 14040-14045 (2005). 

13. Drexler, J. F. et al. Genomic characterization of severe acute respiratory syndrome-related coronavirus in European bats and 
classification of coronaviruses based on partial RNA-dependent RNA polymerase gene sequences. J Virol 84, 11336-11349 (2010). 

14. Lau, S. K. et al. Severe Acute Respiratory Syndrome (SARS) Coronavirus ORF8 Protein Is Acquired from SARS-Related Coronavirus 
from Greater Horseshoe Bats through Recombination. Journal of virology 89, 10532-10547 (2015). 

15. Wu, Z. et al. ORF8-Related Genetic Evidence for Chinese Horseshoe Bats as the Source of Human Severe Acute Respiratory 
Syndrome Coronavirus. The Journal of infectious diseases 213, 579-583 (2016). 

16. Chen, C. Y. et al. Open reading frame 8a of the human severe acute respiratory syndrome coronavirus not only promotes viral 
replication but also induces apoptosis. The Journal of infectious diseases 196, 405-415 (2007). 

17. Wong, H. H. et al. Accessory proteins 8b and 8ab of severe acute respiratory syndrome coronavirus suppress the interferon signaling 
pathway by mediating ubiquitin-dependent rapid degradation of interferon regulatory factor 3. Virology 515, 165-175 (2018). 

18. Keng, C. T. et al. The human severe acute respiratory syndrome coronavirus (SARS-CoV) 8b protein is distinct from its counterpart 
in animal SARS-CoV and down-regulates the expression of the envelope protein in infected cells. Virology 354, 132-142 (2006). 

19. Law, P. Y. et al. Expression and functional characterization of the putative protein 8b of the severe acute respiratory syndrome- 
associated coronavirus. FEBS letters 580, 3643-3648 (2006). 

20. Sung, S. C., Chao, C. Y., Jeng, K. S., Yang, J. Y. & Lai, M. M. The 8ab protein of SARS-CoV is a luminal ER membrane-associated 
protein and induces the activation of ATF6. Virology 387, 402-413 (2009). 





SCIENTIFIC REPORTS | (2018) 8:15177 | DOI:10.1038/s41598-018-33487-8 9 





www.nature.com/scientificreports/ 


21. 


22. 


23. 


24. 


25. 


26. 


27. 


28. 


29. 


30. 
31. 


32. 


33. 
34. 


35. 


36. 


37. 


38. 


39. 


40. 


41. 
42. 


43. 


44. 


45. 


46. 


47. 
48. 


49. 


50. 


51. 


32; 


53. 


54. 


55. 


56. 


57. 


58. 


59. 


60. 


61. 


62. 


Oostra, M., de Haan, C. A. & Rottier, P. J. The 29-nucleotide deletion present in human but not in animal severe acute respiratory 
syndrome coronaviruses disrupts the functional expression of open reading frame 8. Journal of virology 81, 13876-13888 (2007). 
Yount, B. et al. Severe acute respiratory syndrome coronavirus group-specific open reading frames encode nonessential functions 
for replication in cell cultures and mice. Journal of virology 79, 14909-14922 (2005). 

Payne, D. C. et al. Multihospital Outbreak of a Middle East Respiratory Syndrome Coronavirus Deletion Variant, Jordan: A 
Molecular, Serologic, and Epidemiologic Investigation. Open Forum Infect Di 5 (2018). 

Lamers, M. M. et al. Deletion Variants of Middle East Respiratory Syndrome Coronavirus from Humans, Jordan, 2015. Emerging 
infectious diseases 22, 716-719 (2016). 

Pfefferle, S. et al. Reverse genetic characterization of the natural genomic deletion in SARS-Coronavirus strain Frankfurt-1 open 
reading frame 7b reveals an attenuating function of the 7b protein in-vitro and in-vivo. Virology journal 6, 131 (2009). 

Schelle, B., Karl, N., Ludewig, B., Siddell, S. G. & Thiel, V. Selective replication of coronavirus genomes that express nucleocapsid 
protein. Journal of virology 79, 6620-6630 (2005). 

Sola, I., Moreno, J. L., Zuniga, S., Alonso, S. & Enjuanes, L. Role of nucleotides immediately flanking the transcription-regulating 
sequence core in coronavirus subgenomic mRNA synthesis. Journal of virology 79, 2506-2516 (2005). 

McKimm-Breschkin, J. L. & Holmes, I. H. Conditions required for induction of interferon by rotaviruses and for their sensitivity to 
its action. Infect Immun 36, 857-863 (1982). 

Emeny, J. M. & Morgan, M. J. Regulation of the interferon system: evidence that Vero cells have a genetic defect in interferon 
production. The Journal of general virology 43, 247-252 (1979). 

Perlman, S. & Netland, J. Coronaviruses post-SARS: update on replication and pathogenesis. Nat Rev Microbiol 7, 439-450 (2009). 
Biesold, S. E. et al. Type I interferon reaction to viral infection in interferon-competent, immortalized cell lines from the African 
fruit bat Eidolon helvum. PloS one 6, e28131 (2011). 

Kuhl, A. et al. Comparative analysis of Ebola virus glycoprotein interactions with human and bat cells. The Journal of infectious 
diseases 204(Suppl 3), S840-849 (2011). 

Hofmann, A. et al. Epigenetic regulation of lentiviral transgene vectors in a large animal model. Mol Ther 13, 59-66 (2006). 

Ehlen, L. et al. Epithelial cell lines of the cotton rat (Sigmodon hispidus) are highly susceptible in vitro models to zoonotic Bunya-, 
Rhabdo-, and Flaviviruses. Virology journal 13, 74 (2016). 

Lau, S. K. et al. Ecoepidemiology and complete genome comparison of different strains of severe acute respiratory syndrome-related 
Rhinolophus bat coronavirus in China reveal bats as a reservoir for acute, self-limiting infection that allows recombination events. 
Journal of virology 84, 2808-2819 (2010). 

Drosten, C. et al. Identification of a novel coronavirus in patients with severe acute respiratory syndrome. The New England journal 
of medicine 348, 1967-1976 (2003). 

Zwart, M. P. & Elena, S. F. Matters of Size: Genetic Bottlenecks in Virus Infection and Their Potential Impact on Evolution. Annu Rev 
Virol 2, 161-179 (2015). 

Varble, A. et al. Influenza A virus transmission bottlenecks are defined by infection route and recipient host. Cell Host Microbe 16, 
691-700 (2014). 

Bergstrom, C. T., McElhany, P. & Real, L. A. Transmission bottlenecks as determinants of virulence in rapidly evolving pathogens. 
Proceedings of the National Academy of Sciences of the United States of America 96, 5095-5100 (1999). 

Novella, I. S., Elena, S. F., Moya, A., Domingo, E. & Holland, J. J. Size of genetic bottlenecks leading to virus fitness loss is determined 
by mean initial population fitness. Journal of virology 69, 2869-2872 (1995). 

McCrone, J. T. & Lauring, A. S. Genetic bottlenecks in intraspecies virus transmission. Curr Opin Virol 28, 20-25 (2018). 

Xie, Q. et al. Two deletion variants of Middle East respiratory syndrome coronavirus found in a patient with characteristic 
symptoms. Archives of virology 162, 2445-2449 (2017). 

Lu, X. et al. Spike gene deletion quasispecies in serum of patient with acute MERS-CoV infection. Journal of medical virology 89, 
542-545 (2017). 

Niemeyer, D. et al. Middle East respiratory syndrome coronavirus accessory protein 4a is a type I interferon antagonist. J Virol 87, 
12489-12495 (2013). 

Kindler, E. et al. Efficient Replication of the Novel Human Betacoronavirus EMC on Primary Human Epithelium Highlights Its 
Zoonotic Potential. mBio 4 (2013). 

Diehl, W. E. et al. Ebola Virus Glycoprotein with Increased Infectivity Dominated the 2013-2016 Epidemic. Cell 167, 1088-1098 
e1086 (2016). 

Urbanowicz, R. A. et al. Human Adaptation of Ebola Virus during the West African Outbreak. Cell 167, 1079-1087 e1075 (2016). 
Dietzel, E., Schudt, G., Krahling, V., Matrosovich, M. & Becker, S. Functional Characterization of Adaptive Mutations during the 
West African Ebola Virus Outbreak. Journal of virology 91 (2017). 

Ueda, M. T. et al. Functional mutations in spike glycoprotein of Zaire ebolavirus associated with an increase in infection efficiency. 
Genes Cells 22, 148-159 (2017). 

Wang, M. K., Lim, S. Y., Lee, S. M. & Cunningham, J. M. Biochemical Basis for Increased Activity of Ebola Glycoprotein in the 
2013-16 Epidemic. Cell Host Microbe 21, 367-375 (2017). 

Marzi, A. et al. Recently Identified Mutations in the Ebola Virus-Makona Genome Do Not Alter Pathogenicity in Animal Models. 
Cell Rep 23, 1806-1816 (2018). 

Balboni, A., Palladini, A., Bogliani, G. & Battilani, M. Detection of a virus related to betacoronaviruses in Italian greater horseshoe 
bats. Epidemiol Infect 139, 216-219 (2011). 

Rihtaric, D., Hostnik, P., Steyer, A., Grom, J. & Toplak, I. Identification of SARS-like coronaviruses in horseshoe bats (Rhinolophus 
hipposideros) in Slovenia. Archives of virology 155, 507-514 (2010). 

Vallo, P., Guillen-Servent, A., Benda, P., Pires, D. B. & Koubek, P. Variation of mitochondrial DNA in the Hipposideros caffer 
complex (Chiroptera: Hipposideridae) and its taxonomic implications. Acta Chiropterol 10, 193-206 (2008). 

Eckerle, I. et al. Replicative Capacity of MERS Coronavirus in Livestock Cell Lines. Emerging infectious diseases 20, 276-279 (2014). 
Eckerle, I. et al. Bat airway epithelial cells: a novel tool for the study of zoonotic viruses. PloS one 9, e84679 (2014). 

Dijkman, R. et al. Isolation and characterization of current human coronavirus strains in primary human epithelia cultures reveals 
differences in target cell tropism. Journal of virology (2013). 

Muller, M. A. et al. Human coronavirus EMC does not require the SARS-coronavirus receptor and maintains broad replicative 
capability in mammalian cell lines. MBio 3 (2012). 

Muller, M. A. et al. Human coronavirus NL63 open reading frame 3 encodes a virion-incorporated N-glycosylated membrane 
protein. Virology journal 7, 6 (2010). 

Hofmann, H. et al. Susceptibility to SARS coronavirus S protein-driven infection correlates with expression of angiotensin 
converting enzyme 2 and infection can be blocked by soluble receptor. Biochem Biophys Res Commun 319, 1216-1221 (2004). 
Morris, D. W. et al. An inherited duplication at the gene p21 Protein- Activated Kinase 7 (PAK7) is a risk factor for psychosis. Hum 
Mol Genet 23, 3316-3326 (2014). 

Matrosovich, M., Matrosovich, T., Garten, W. & Klenk, H. D. New low-viscosity overlay medium for viral plaque assays. Virology 
journal 3, 63 (2006). 





SCIENTIFIC REPORTS | (2018) 8:15177 | DOI:10.1038/s41598-018-33487-8 10 


www.nature.com/scientificreports/ 


Acknowledgements 

We would like to thank Artem Siemens, Stephan Kallies, Tobias Bleicker, Sebastian Briinink, and Monika 
Eschbach-Bludau for excellent technical assistance. We are grateful to Matthias Lenk (Friedrich-Loeffler-Institut, 
Greifswald-Insel Riems) for providing goat and sheep cell lines. This study was supported by the Deutsche 
Forschungsgemeinschaft (DFG DR 772/12-1) and EU-funded project COMPARE GA no. 643476. Furthermore, 
we acknowledge support from the Deutsche Forschungsgemeinschaft (DFG) and the Open Access Publication 
Fund of Charité - Universitatsmedizin Berlin. 


Author Contributions 
C.D., M.A.M., D.M. conceived and designed the experiments. V.M.C., EG.R., A.B., M.B., D.R., I.T., R.S.A. and 
J.F.D. sampled and processed bat specimens. T.B., H.R. and L.T.G. established and characterized transgenic 
bat cell line. R.D., V.T. and A.P. provided essential materials. D.M. and L.T.G. performed cloning and infection 
experiments with recombinant viruses. D.M., C.D. and M.A.M. interpreted the data and wrote the manuscript 
with input from all authors. 


Additional Information 
Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-33487-8. 


Competing Interests: The authors declare no competing interests. 


Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and 
institutional affiliations. 


Open Access This article is licensed under a Creative Commons Attribution 4.0 International 


License, which permits use, sharing, adaptation, distribution and reproduction in any medium or 
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- 
ative Commons license, and indicate if changes were made. The images or other third party material in this 
article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the 
material. If material is not included in the article’s Creative Commons license and your intended use is not per- 
mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the 
copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. 


© The Author(s) 2018 





SCIENTIFIC REPORTS | (2018) 8:15177 | DOI:10.1038/s41598-018-33487-8 11 


